Multimodal transition and excitability of a neural oscillator 
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We analyze the response of the Morris-Lecar model to a periodic train of short current pulses in the 
period-amplitude plane. For a wide parameter range encompassing both class 2 and class 3 behavior 
in Hodgkin's classification there is a multimodal transition between the set of odd modes and the 
set of all modes. It is located between the 2:1 and 3:1 locked- in regions. It is the same dynamic 
instability as the one discovered earlier in the Hodgkin-Huxley model and observed experimentally 
in squid giant axons. It appears simultaneously with the bistability of the states 2:1 and 3:1 in 
the perithreshold regime. These results imply that the multimodal transition may be a universal 
property of resonant neurons. 
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I. INTRODUCTION 

In 1948 Hodgkin studied the response of neurons to 
stimulation by a constant current [1]. He summarized 
his findings by dividing neurons into three classes: those 
having continuous relation between the current ampli- 
tude and response frequency (type 1), those with a dis- 
continuous jump of response frequency at the stimulus 
threshold (type 2), and those which spike only once or 
twice to the constant current stimulus (type 3). It is 
now a widely accepted view that integrator (resonator) 
neurons exhibit type 1 (type 2) behavior, respectively. 
Well known examples of type 2 neurons include Hodgkin- 
Huxley (HH) [2J, fast-spiking cortical cells [sl, 0], Morris- 
Lecar (ML) [5|,l6| and Hindmarsh-Rose models. Some 
neuron models are known to exhibit different types of 
excitability, depending on parameter values. Prescott et 
al. [H showed in the ML model, originally used to de- 
scribe the barnacle giant muscle fiber, that change of one 
parameter was sufficient to switch between type 1, type 
2, and type 3 dynamics. Mathematically speaking, class 
1 excitability occurs through a saddle-node bifurcation, 
class 2 excitability occurs through a Hopf bifurcation, 
and class 3 excitability is a result of a quasiseparatrix 
crossing. The ML model consists of equations describing 
the time evolution of two dynamic variables, the fast and 
the slow one. Prescott et al. [s] emphasized that the 
change of excitability in this system was a result of com- 
petition between slow and fast currents. The ML model 
was also used as a minimal model in the investigation of 
changes of the dynamics of the CAl pyramidal neurons 
between type 1 and type 2 [9]. 

There is more evidence that the spike initiation mech- 
anism is not a fixed property of the neuron. In some 
experiments the squid giant axons had type 3 instead of 
type 2 excitability [Tol, |Tl| . The discrepancy between the 
type 2 behavior of the HH model and experiment was 
explained by modifying a single parameter in the term 
describing the potassium current In a recent study 
of the periodically stimulated HH model by the present 
author it was found that the firing rate may be either con- 
tinuous or discontinuous function of the current ampli- 



tude, depending on the stimulus frequency [Ts*]. Bistable 
behavior at the excitation threshold was identified as a 
signature of an antiresonant regime of a resonant neuron. 
When brief stimuli arrive at resonant frequencies, the HH 
neuron may respond with arbitrarily low firing rate. The 
dependence of the firing rate on the current amplitude 
scales with a square root above the threshold, consistent 
with a saddle- node bifurcation p^ . 

The HH neuron's response at resonant frequencies can 
be divided into three regimes: (i) short pulses, where the 
width r does not exceed the optimal width Topt associ- 
ated with a minimum threshold Topt^ (ii) Topt < r < Tresi 
and r Tres [13]. The Hodgkin classification scheme 
is related to case (iii). The analysis of response in the 
limit of short stimuli gives a much more interesting in- 
formation about the neuron's dynamics [13]. Since the 
behavior of the HH model at the threshold depends on 
stimulation frequency, it would be useful to investigate 
the same problem in a simpler neural oscillator, such as 
the ML model, which can be viewed as a minimal model 
of a resonant neuron. 

An important question related to the study of ex- 
citability is the existence of the dynamic singularity in 
the ML model. The recently discovered multimodal tran- 
sition (MMT) [l5| in the response of the HH cells at fre- 
quencies above the main resonant frequency manifests it- 
self as a change of the parity of response modes, strongly 
affecting the dynamics in its vicinity. Experimental data 
of Takahashi et al. [16] provide strong evidence for the 
existence of this transition [l^. The MMT occurs just 
above the threshold, between the locked- in states 2:1 and 
3:1. Could this be a universal property of resonant neu- 
rons? How does the MMT relate to the Hodgkin's clas- 
sification? 

The knowledge of the perithreshold behavior of var- 
ious groups of neurons is crucial if we want to under- 
stand their role in encoding different types of neural in- 
put [sl, 3, [m, . There is also another important reason 
for studying patterns of excitability. Finding regimes of 
similar behavior in different models would provide jus- 
tification for replacing more complex sets of equations 
by simpler ones. For example, we can ask under what 
circumstances the ML model, or its equivalent, could be 
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used as a reasonable simplification of the HH-type mod- 
els. It was shown earlier that the relaxation times of both 
the HH model and the type 2 version of the ML model 
just below the threshold of repetitive firing, when stimu- 
lated by a constant current, were both equal to 1/2 [20]. 
Does this similarity extend to finite stimulus frequencies? 

II. THE MODEL AND RESULTS 

We use the form of the ML model proposed by Prescott 
et al. [8], 
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FIG. 1: Response diagram for — 0. The main locked- in 
state are labeled by the inverse of the firing rate. 
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The fast activation variable V competes with the slow re- 
covery variable w. Parameter values were chosen in Ref. 
m to produce different spiking patterns: ENa = 50mV, 
Ek -lOOmV, El -70mV, gfast = 20mS/cm^ 
Qsiow = 20mS/cm^, gL 2mS/cm^, (/) = 0.15, f3m = 
-1.2mV, jm = 18mV, and 7^ = lOmV. C 2/iF/cm^ 
is the membrane capacitance. We chose the input cur- 
rent to be a periodic set of rectangular steps of period 
Ti, height Iq and width r = 0.5ms. Studying the HH 
model, we learned that the topology of the global bifur- 
cation diagram is only weakly dependent on shape details 
of individual pulses, provided they remain short com- 
pared to the time scale of the main resonance [13|, 121} . 
The calculations are carried out within the fourth-order 
Runge-Kutta scheme with the time step of 0.001ms. In- 
dividual runs at fixed parameters were carried out for 
lOOOT^. Since the variation of is sufficient to alter the 
excitability type of the model, we study the effect of /3u) 
on the dynamics at finite frequencies. Changes of other 
parameters, Prn, gfast, gsiow, 7m, and 7^, may result in 
similar evolution of the neuron's dynamics [8]. 

Figure [T] shows the response diagram in the period- 
amplitude plane for = 0. In the limit /q = const 
this choice leads to type 1 excitability [8|. In Fig. [1] 



the dependence of the firing rate fo/fi on Iq is contin- 
uous everywhere along the excitation threshold, scaling 
approximately as {Iq — IthY^'^, where Ith is the value of 
Iq at the threshold. 








FIG. 2: Response diagram for = — 13mV. The locked-in 
states of order 3:1 and higher are bistable along the threshold. 
However the 2:1 state remains monostable. The dashed line 
separates the monostable regime from the bistable area. We 
have verified that MMT does not appear in this case. 

When = — 13mV, the neuron displays class 2 dy- 
namics for a constant current. In Fig. [2] we can see 
that the firing rate is a discontinuous function of the cur- 
rent amplitude also at high stimulation frequencies. For 
Ti < 4ms almost the entire threshold is bistable. How- 
ever, for Ti > 4ms the system remains monostable. More 
precisely, the bistability does not extend beyond the 3:1 
state and the edge of the 2:1 state is monostable. The 
low- frequency j"esponse in Figs. [T] and [2] is very similar. 
Wang et al. [22| also noted the similarity of response in 
this regime for a sinusoidal input. 

Fig. [3] shows the response diagram for = — 23mV. 
There are now large bistable regions of the states 3:1 
and 2:1. The high-frequency part of the diagram, for 
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FIG. 3: Response diagram for I3w = — 23mV. The dashed 
hne separates the monostable regime from the bistable area. 
The location of the odd-all multimodal transition is denoted 
by full squares. For 1.5ms < Ti < 5ms this diagram closely 
resembles the one obtained for the HH model 13]. 
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FIG. 4: The firing rate (top) and the weight of even modes 
(bottom) for fiw — — 23mV and Iq — 245/iA/cm^. 



Ti < 5ms, closely resembles the regime of the HH model 
where the MMT occurs [H, [I^ • We have analyzed the 
ISI histogram and found the same dynamic singularity 
in the ML model. The location of the transition is indi- 
cated in Fig. [3] by full squares. The fo vs Iq dependence 
in the interval between the MMT and the 2:1 state is ap- 
proximately linear, as in the HH model [15]. The MMT 
along the Ti axis for = — 23mV is shown in Fig. IH 
The weight of even modes drops sharply in the vicinity of 
the minimum of the firing rate. The edges of individual 
modes near the transition scale logarithmically, as in the 
HH model [iBi] . All these signatures of the MMT and the 
topology of the bifurcation diagram in the vicinity of the 
MMT in the ML model are identical to the HH model. 
Therefore it is reasonable to conclude that (i) these may 
be universal properties associated with the MMT, and 
(ii) MMT occurs in many other type 2 or type 3 reso- 
nant neurons. 

Fig. [5] shows sample V{t) runs on both sides of 
the even-all transition. In the top panel, obtained at 
Ti = 2.65ms, the even modes 4:1 and 8:1 dominate. In 
the bottom panel obtained at Ti = 2.45ms only odd 
modes are present. Here the height of the V{t) peak 
correlates with the length of the preceding ISI. A careful 
examination of the peak heights reveals a preference for 
a To = 3T^, following the action potential with the max- 
imum value of V. This preference is manifested either as 
(i) another action potential, or as (ii) subthreshold peak 
with a height somewhat larger than that of its immediate 
neighbors. Judging by the height of subthreshold peaks 
there is a clear preference for a period 2T^ subthresh- 
old oscillation. We can view the odd-only periods in the 
bottom panel as a sum of 3T^ + 2nTi, where n = 0, 1, — 

We also analyzed the ISI histogram as a function of 
Z^^y, looking for signs of the MMT. For = — 13mV we 
can see a precursor of the odd-all MMT in Fig. [6l There 
is a local minimum of fo/ fi at Iq = 198.2/iA/cm^ and 
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FIG. 5: Sample V{t) run for = -23mV, lo = 245/iA/cm^. 
The top panel contains data taken above the odd-all transition 
at T = 2.65ms. Here the even modes are more frequent than 
the odd ones but there is also some admixture of the odd 
modes. The bottom panel shows data at Ti = 2.45ms. Note 
the absence of even multiples of the driving period in the 
bottom panel. Here the response contains only the following 
multiples of Ti: 3, 5, 7, and 9. 



a significant decrease of the participation rate of even 
modes close to the threshold. 

The MMT is tied to the appearance of bistability 
along the bottom edge of the 2:1 state, which occurs 
for < — 14.5mV. In Fig. [7] we can see that for 
/3yj = —18ms the MMT is already well developed. As 
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FIG. 6: The firing rate and tiie joint histogram weight of even 
modes for I3w = — 13mV and Ti = 4.7ms. A precursor of the 
multimodal transition can be seen near Iq = 198.2/iA/cm^. 
There is a well pronounced local minimum of the fo/ fi and a 
significant reduction of the participation rate of even modes. 



the MMT is approached from above, only very high or- 
der even modes remain. The edges of the even modes 
scale logarithmically as a function of |/o — ^*|, where 
/* is the current amplitude at which the transition oc- 
curs. When the bist ability of the state 2:1 disappears 
near = — 29mV, so does the MMT. At this value of 
/3w the 3:1 state vanishes from the global bifurcation di- 
agram. 



III. CONCLUSIONS 

The dynamics of neurons at frequencies above the res- 
onant frequency is not directly related to either class 2 
or class 3 excitability. The global bifurcation diagram of 
the class 3 ML model for (3^^ = — 23mV closely resembles 
the class 2 HH model [l5| at intermediate frequencies. 
This is consistent with the remarks of Prescott et al. [H 
that neurons should not be labeled as being strictly type 
2 or type 3. The transition to bistability at the edge of 
a 2:1 state in a model stimulated by a train of current 
pulses does not occur for the same parameter values as 
the transition to bistability for a constant current. The 
bistability appears first in the limit of high frequency 
stimulation. As P^j decreases further, the perithreshold 
regions for larger Ti also become bistable. The MMT 
occurs when both the 2:1 and the 3:1 state have bistable 
regions along the threshold. The ability to predict the 
existence of the MMT from the topology of the locked-in 
states in the perithreshold regime is of course a very use- 
ful property. Testing for bistabilities near the threshold 
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FIG. 7: From top to bottom: the firing rate, the joint 
histogram weight of even modes, and the ISI spectrum for 
— — 18mV and Ti — 3.7ms. The multimodal transition 
occurs at Iq — 207.5/iA/cm^. 



is much simpler and computationally much more efficient 
than analyzing the evolution of the entire ISI histogram 
as a function of pulse frequency or amplitude. 

Wang et al. ^] noted the similarity of class 1 and 
class 2 neuron response to sinusoidal signal at low fre- 
quencies. This is not surprising, once we note, that the 
threshold bistability, and a discontinuity of the fo/fi 
associated with it, appears first in the limit of high fi 
and extends towards low fi with decreasing P^. When 
< — 14.5mV, the 3:1 state becomes bistable and the 
entire threshold at high frequencies rises significantly. 
The emergence of class 3 behavior may be viewed as a 
strong rigidity of neuron dynamics at high fi. 

We could also think of an alternative classification of 
neuron excitability for stimuli of finite frequency: (i) class 
1, where the firing rate is a continuous function of /q 
everywhere along the threshold, (ii) class 2, with bista- 
bilities at the edges of high-order locked-in states, (iii) 
class 3, where the MMT exists, and (iv) class 4, where 
both the MMT and the bistabilities are absent and the 
neuron responds to a constant current by emitting only 
a few spikes. 

The presence of the MMT in both the ML and HH 
model suggests that the same transition is present in 
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other resonant neuron models classified as type 2 or type in Gdansk. 
3 cells. The existence of MMT may have important phys- 
iological consequences. It may be relevant in the high- 
frequency auditory nerve fiber stimulation [23] and pos- 
sibly in the clinical procedure of deep brain stimulation 
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